Stable Computation of Pad e - Hermite and Simultaneous Pad
نویسندگان
چکیده
In a recent paper Cabay et al. 1996a], the authors develop a fast, iterative, look-ahead algorithm for numerically computingPad e-Hermite systems and simultaneous Pad e systems along a diagonal of the associatedPad e tables. Included there is a detailed error analysis showing that the algorithm is weakly stable. In this paper, we describe a Fortran implementation, VECTOR PADE, of this algorithm together with a number of numerical experiments. These experiments show that the theoretical error bounds obtained in Cabay et al. 1996a] reeect the general behavior of the actual error, but that in practice these bounds are large overestimates. Permission to make digital or hard copies of part or all of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for proot or direct commercial advantage and that copies show this notice on the rst page or initial screen of a display along with the full citation. Copyrights for components of this work owned by others than ACM must be honored. Abstracting with credit is permitted. To copy otherwise, to republish, to post on servers, to redistribute to lists, or to use any component of this work in other works, requires prior speciic permission and/or a fee. Permissions may be requested from Publications
منابع مشابه
Recursiveness in Matrix Rational Interpolation Problems 1 2
We consider the problem of computing solutions to a variety of matrix rational interpolation problems. These include the partial realization problem for matrix power series and Newton-Pad e, Hermite-Pad e, Simultaneous Pad e, M-Pad e and multipoint Pad e approximation problems along with their matrix generalizations. A general recurrence relation is given for solving these problems. Unlike othe...
متن کاملInversion Components of Block Hankel-like Matrices
The inversion problem for square matrices having the structure of a block Hankel-like matrix is studied. Examples of such matrices in&de Hankel striped, Hankel layered, and vector Hankel matrices. It is shown that the components that both determine nonsingularity and construct the inverse of such matrices are closely related to certain matrix polynomials. These matrix polynomials are multidimen...
متن کاملStrong asymptotics of orthogonal polynomials with varying measures and Hermite–Pad e approximants
The strong asymptotic behaviour of orthogonal polynomials with respect to a general class of varying measures is given for the case of the unit circle and the real line. These results are used to obtain certain asymptotic relations for the polynomials involved in the construction of Hermite–Pad e approximants of a Nikishin system of functions. c © 1998 Elsevier Science B.V. All rights reserved....
متن کاملMultiple orthogonal polynomials
Results on multiple orthogonal polynomials will be surveyed. Multiple orthogonal polynomials are intimately related to Hermite–Pad e approximants and often they are also called Hermite–Pad e polynomials. Special attention will be paid to an application of multiple orthogonal polynomials and to analytic theory of two model families of general multiple orthogonal polynomials, referred to as Angel...
متن کامل